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I The quark distributions in the proton and neutron are parametrized as Parton Distri- 

^ ■ bution Functions (PDFs) obtained from global fits to various sets of data at very high 

O ■ energies. These fits are done within the theory of Quantum Chromodynamics (QCD) in 

^ I either leading order (LO) or next to leading order (NLO). The most important data come 

from deep-inelastic e//i scattering experiments on hydrogen and deuterium, and and 
^ . experiments on nuclear targets. In previous publications d 0,11] we have compared 

I I the predictions of the NLO MRSR2 PDFs to deep-inelastic e//i scattering data [4] on 

hydrogen and deuterium from SLAC, BCDMS and NMC. In order to get agreement 
^ . with the lower energy SLAC data for F2 and R down to Q^=l GeV^, and at the highest 

values of X (x = 0.9), we found that the following modifications to the NLO MRSR2 
^ PDFs must be included. 

. 1 . The relative normalizations between the various data sets and the BCDMS system 

atic error shift must be included M 



2. Deuteron binding corrections need to be applied and the ratio of d/u at high x must 
be increased as discussed in ref. [1]. 

3. Kinematic higher- twist originating from target mass effects |0] are very large and 
must be included. 

4. Dynamical higher-twist corrections are smaller but also need to be included 

5. In addition, our analysis including QCD Next to NLO (NNLO) terms shows 10] 
that most of the dynamical higher-twist corrections needed to fit the data within a 
NLO QCD analysis originate from the missing NNLO higher order terms. 

Our analysis shows that the NLO MRSR2 PDFs with target mass and NNLO higher 
order terms describe electron and muon scattering F2 and R data with a very small 
contribution from higher twists. Studies by other authors [6] also show that in NNLO 
analyses the dynamic higher twist corrections are very small. If (for > I GeV^) most 



of the higher-twist terms needed to obtain agreement with the low energy data actually 
originate from target mass effects and missing NNLO terms (i.e. not from interactions 
with spectator quarks) then these terms should be the same in and e//i scattering. 
Therefore, low energy data should be described by the PDFs which are fit to high 
energy data and are modified to include target mass and higher-twist corrections that fit 
low energy djl scattering data. However, for < I GeV^ additional non-perturbative 
effects from spectator quarks must also be included \J^. 

PREVIOUS RESULTS WITH GRV94 PDFS AND Xw 

In a previous communication ^ we used a modified scaling variable Xw and fit for 
modifications to the GRV94 leading order PDFs such that the PDFs describe both 
high energy low energy e//i data. In order to describe low energy data down to the 
photoproduction limit (Q^ = 0), and account for both target mass and higher twist 
effects, the following modifications of the GRV94 LO PDFs are need: 

1. We increased the d/u ratio at high x as described in our previous analysis Ol]. 

2. Instead of the scaling variable x we used the scaling variable x^ = (Q^ + 
B)/{2Mv +A) (or =x{Q^ +B)/{Q^ +Ax)). This modification was used in early fits 
to SLAC data |9]. The parameter A provides for an approximate way to include 
both target mass and higher twist effects at high x, and the parameter B allows the 
fit to be used all the way down to the photoproduction limit (2^=0). 

3. In addition as was done in earlier non-QCD based fits [ 10] to low energy data, we 
multiplied all PDFs by a factor K=Q^ I {Q^ -\-C). This was done in order for the 
fits to describe low data in the photoproduction limit, where F2 is related to the 
photoproduction cross section according to 

^(ri^) = q2 ^2 = ^2 ^2 

4. Finally, we froze the evolution of the GRV94 PDFs at a value of = 0.24 (for 

< 0.24), because GRV94 PDFs are only valid down to = 0.23 GeV^. 

In our analyses, the measured structure functions were corrected for the BCDMS 
systematic error shift and for the relative normalizations between the SLAC, BCDMS 
and NMC data llliS]. The deuterium data were corrected for nuclear binding effects lH 
0]. A simultaneous fit to both proton and deuteron SLAC, NMC and BCDMS data (with 
X > 0.07 only) yields A=1.735, B=0.624 and C=0.188 (GeV^) with GRV94 LO PDFs 
(X^ = 1351/958 DOF). Note that for the parameter A accounts for both target mass 
and higher twist effects. 



NEW ANALYSIS WITH ^w, Go AND GRV98 PDFS 

In this publication we update our previous studies, which were done with a new 
improved scaling variable and fit for modifications to the more modem GRV98 LO 




FIGURE 1. Electron and muon F2 data (SLAC, BCDMS, NMC, HI 94) used in our GRV98 ^v, fit 
compared to the predictions of the unmodified GRV98 PDFs (LO, dashed Hne) and the modified GRV98 
PDFs fits (LO+HT, solid line); [a] for F2 proton, [b] for F2 deuteron, and [c] for the HI and NMC proton 
data at low x. 



PDFs such that the PDFs describe both high energy and low energy electron/muon data. 
We now also include NMC and HI 94 data at lower x. Here we freeze the evolution 
of the GRV98 PDFs at a value of = (for < 0.8), because GRV98 PDFs 
are only valid down to =0.8 GeV^. In addition, we use different photoproduction 
limit multiplicative factors for valence and sea. Our proposed new scaling variable 
is based on the following derivation. Using energy momentum conservation, it can 
be shown that the factional momentum ^ = (p^ + po) / {Pz + Po) carried by a quark 
of 4-mometum p in a proton target of mass M and 4-momentum P is given hy ^ = 
xQ'^/[0.5Q^{\ + [1 + {2Mx)^/Q^]^/% where 

22'^ = [2^ + M/ -M,2] + [(^2 _^M/ -M,2)2 + 4e2(M,2 + /'2)]i/2. 

Here Mj is the initial quark mass with average initial transverse momentum Pj and 




FIGURE 2. Comparisons to proton and iron data not included in our GRV98 fit. (a) Comparison 
of SLAC and JLab (electron) F2p data in the resonance region (or fits to these data) and the predictions 
of the GRV98 PDFs with (LO+HT, solid) and without (LO, dashed) our modifications, (b) Comparison 
of photoproduction data on protons to predictions using our modified GRV98 PDFs. (c) Comparison of 
representative CCFR and charged-current differential cross sections 01 [3 on iron at 55 GeV and 
the predictions of the GRV98 PDFs with (LO+HT, solid) and without (LO, dashed) our modifications. 



Mf is the mass of the quark in the final state. The above expression for was previously 
derived Jsl] for the case of Pt = 0. Assuming M,- = we use instead: 
= x{Q^+B+M/)/{0.5Q^{l + [1 + {2Mx)^/Q^] i/^) +Ax) 

Here Mf=0, except for charm-production processes in neutrino scattering for which 
Mf=l.5 GeV. For the parameter A is expected to be much smaller than for since 
now it only accounts for the higher order (dynamic higher twist) QCD terms in the form 
of an enhanced target mass term (the effects of the proton target mass are already taken 
into account using the exact form in the denominator of ^w)- The parameter B accounts 
for the initial state quark transverse momentum and final state quark effective AM 
(originating from multi-gluon emission by quarks). 

Using closure considerations |11] (e.^.the Gottfried sum rule) it can be shown that, 
at low Q^, the scaling prediction for the valence quark part of F2 should be multiplied 
by the factor K=[1-GI(Q^)][1+M(Q^)] where Go = 1/(1+2^/0.71)2 is the proton elastic 
form factor, and M{Q^) is related to the magnetic elastic form factors of the proton and 
neutron. At low Q^, [l-Gl(Q^)] is approximately Q^/(Q^ +C) with C = 0.71/4 = 0.178 
(versus our fit value C=0.18 with GRV94). In order to satisfy the Adler Sum rule lfl2l] 
we add the function M{Q^) to account for terms from the magnetic and axial elas- 
tic form factors of the nucleon). Therefore, we try a more general form ^vaZe;3ce=[l- 
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FIGURE 3. Comparisons to data on deutrerium which were not included in our GRV98 fit. (a) 
Comparison of SLAC and JLab (electron) F2ci data in the resonance region and the predictions of the 
GRV98 PDFs with (LO+HT, solid) and without (LO, dashed) our modifications, (b) Comparison of 
photoproduction data on deuterium to predictions using our modified GRV98 PDFs (including shadowing 
corrections), (c) The shadowing corrections that were applied to the PDFs for predicting the photoproduc- 
tion cross section on deuterium. 



Gl{Q^)][Q^+C2MQ^ +Civ], and K,,a=Q^KQ^+Csea). Using this form with the GRV98 
PDFs (and now also including the very low x NMC and HI 94 data in the fit) we 
find A=0.419, 5=0.223, and Cu,=0.544, C2v=0.431, and C,e„=0.380 (all in GeV^, = 
1235/1200 DOF). As expected, A and B are now smaller with respect to our previous 
fits with GRV94 and Xw With these modifications, the GRV98 PDFs must also be mul- 
tiplied by A^=1.01 1 to normalize to the SLAC data. The fit (Figure[T]) yields the fol- 
lowing normalizations relative to the SLAC F2p data (5LACz)=0.986, BCDMSp=Q.96A, 
BCDMSd=0.9M, NMCp=l. 00, NMCd=0.993, H\p=0.911, and BCDMS systematic er- 
ror shift of 1 .7). (Note, since the GRV98 PDFs do not include the charm sea, for Q^>0.% 
GeV^ we add charm production using the photon-gluon fusion model in order to fit the 
very high v HERA data. This is not needed for any of the low energy comparisons but 
is only needed to describe the highest v HERA electro and photoproduction data). 

Comparisons of predictions using these modified GRV98 PDFs to other data which 
were not included in the fit is shown in Figures El and El From duality ifljll consider- 
ations, with the (^vi scaling variable, the modified GRV98 PDFs should also provide a 
reasonable description of the average value of F2 in the resonance region. Figures 12 a) 
and |3ta) show a comparison between resonance data (from SLAC and Jefferson Lab, 
or parametrizations of these data lIlSll ) on protons and deuterons versus the predictions 
with the standard GRV98 PDFs (LO) and with our modified GRV98 PDFs (LO-i-HT). 
The modified GRVB98 PDFs are in good agreement with SLAC and JLab resonance 
data down to = 0.07 (although resonance data were not included in our fits). There 



is also very good agreement of the predictions of our modified GRV98 in the =0 
limit with photoproduction data on protons and deuterons as shown in Figure |2tb) and 
Etb). In predicting the photoproduction cross sections on deuterium, we have applied 
shadowing corrections [19] as shown in Figure |3tc). We also compare the predictions 
with our modified GRV98 PDFs (LO+HT) to a few representative high energy CCFR 
and charged-current differential cross sections [3, 13] on iron (neutrino data were 
not included in our fit). In this comparison we use the PDFs to obtain F2 and xFt, and 
correct for nuclear effects in iron The structure function 2xFi is obtained by using 
the Rworid fit from reference f?]. There is very good agreement of our predictions with 
these neutrino data on iron. 

In order to have a full description of^ charged current and processes, the 
contribution from quasielastic scattering 1 16] must be added separately at x = 1 . The best 
prescription is to use our model in the region above the first resonance (above W=135 
GeV) and add the contributions from quasielastic and first resonance [17] (VK=1 .23 GeV) 
separately. This is because the W = M and W=l.23 GeV regions are dominated by 
one and two isospin states, and the arnplitudes for neutrino versus electron scattering 
are related via Clebsch-Gordon rules [17] instead of quark charges (also the V and A 
couplings are not equal at low W and Q^). In the region of higher mass resonances 
(e.g. W=l.l GeV) there is a significant contribution from the deep-inelastic continuum 
which is not well modeled by the existing fits iIitIi to neutrino resonance data (and using 
our modified PDFs should be better). For nuclear targets, nuclear corrections |0] must 
also be applied. Recent results from Jlab indicate that the Fe/D ratio in the resonance 
region is the same as the Fe/D ratio from DIS data for the same value of (or t,yj). 
The effects of terms proportional to the muon mass and F4 and^ structure functions 
in neutrino scattering are small and are discussed in Ref. [T^, 18]. In the future, we 
plan to investigate the effects of including the initial state quark Pj in and institute 
further improvements such as allowing for different higher twist parameters for u, d, s, 
c, b quarks in the sea, and the small difference (expected in the Adler sum rule) in the K 
factors for axial and vector terms in neutrino scattering. In addition, we can multiply the 
PDFs by a modulating function A(W,2^) to improve modeling in the resonance 

region (for hydrogen) by including (instead of predicting) the resonance data [ 15] in the 
fit. We can also include resonance data on deuterium lllSh and heavier nuclear targets in 
the fit, and low energy neutrino data. Note that because of the effects of experimental 
resolution and Fermi motion [ 20] (for nuclear targets), a description of the average cross 
section in the resonance region is sufficient for most neutrino experiments. 



In leading order QCD (e.g. GRV98 LO PDFs), F2 for the scattering of electrons and 
muons on proton (or neutron) targets is given by the sum of quark and anti-quark 
distributions (each weighted the square of the quark charges): 



APPENDIX 



F2{x) 
2xFi (jc) 



i:iej[xqi{x)+xqi{x)] 
F2{x){l+4MxyQ^)/{l+K,rld)- 



(1) 
(2) 



Here, R^orid is parameterized 10] by: 

/ 0-0635 ^, 0.5747 0.3534 

i?woW.(x,e) - /^^(22/o.04)^^^'2^ + ^ 24 + 0.09' 

where = l. + ^x^^^. 

The Rworid function provides a good description of the world's data in the > 
0.5 and x > 0.05 region. Note that the R^orid function breaks down below 2^ = 0.3. 
Therefore, we freeze the function at = 0.35 and introduce the following function 
for R in the < 0.35 region. The new function provides a smooth transition from 
= 0.35 down to 2^ = by forcing R to approach zero at 2^ = as expected in 
the photoproduction limit (while keeping a 1 /2^ behavior at large and matching to 
Rworid at 2^ = 0.35). 

R{x,Q^) = 3.207 x^^x7?wrf(x,2' = 0.35). (4) 

In the comparison with CCFR charged-current differential cross section on iron, 
a nuclear correction for iron targets is applied. We use the following parameterized 
function, f{x) (fit to experimental electron and muon scattering data for the ratio of 
iron to deuterium cross sections), to convert deuterium structure functions to (isoscalar) 
iron structure functions 

f{x) - 1.096 -0.364;c-0.278e"2i.94x + 2.772;ci4.4i7 ^^^^ 

For the ratio of deuterium cross sections to cross sections on free nucleons we use 
the following function obtained from a fit to SLAC data on the nuclear dependence of 
electron scattering cross sections |3]. 

f{x) = (0.985 ±0.0013) X (1+0.422jc-2.745jc2 

+7.570;c3- 10.335/ + 5.422^5). (6) 

This correction is only valid in the 0.05 < x < 0.75 region. In neutrino scattering, we 
use the same nuclear correction factor for F2, XF2, and 2xF\ . 

The d/u correction for the GRV98 LO PDFs is obtained from the NMC data for 
F2 ■ Here, Eq. |6lis used to remove nuclear binding effects in the NMC deuterium 
F2 data. The correction term, d{d/u) is obtained by keeping the total valence and sea 
quarks the same. 

5{d/u){x) = -0.008 17 + 0.0506jc + 0.0798jc2, (7) 



where the corrected d/u ratio is {d/u)' = (d/u) + 5{d/u). Thus, the modified u and d 
valence distributions are given by 



The same formalism is applied to the modified u and d sea distributions. Accidently, the 
modified u and d sea distributions (based on NMC data) agree with the NUSEA data 
in the range of x between 0. 1 and 0.4. Thus, we find that any futher correction on sea 
quarks is not necessary. 
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